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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si)  AND 
METRIC  (SI)  TO  U.  S.  CUSTOMARY  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can 
follows : 


Multiply By 

U.  S.  Customary  to  Metric  (SI 


feet 

feet  per  second  per 
second 

miles  per  hour  (U.  S. 
statute) 

tons  (2,000  lb,  mass) 


0.301*8 

0.30U8 

1.6093UU 

907.161*7 


Metric  (SI)  to  U. 


S . Customary 


metres 


3.280839 


be  converted  as 

To  Obtain 

l 

metres 

metres  per  second  per 
second 

kilometres  per  hour 
kilograms 

feet 
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A NUMERICAL  MODEL  FOR  TSUNAMI  INUNDATION 


PART  I:  INTRODUCTION 


Tsunamis 


1.  Of  all  water  waves  that  occur  in  nature,  one  of  the  most 
destructive  is  the  tsunami.  The  term  "tsunami,"  originating  from  the 
Japanese  words  "tsu"  (harbor)  and  "nami"  (wave),  is  used  to  describe  sea 
waves  of  seismic  origin.  Tectonic  earthquakes,  i.  e.  earthquakes  that 
cause  a deformation  of  the  seabed,  appear  to  be  the  principal  seismic 
mechanism  responsible  for  the  generation  of  tsunamis.  Coastal  and  sub- 
marine landslides  and  volcanic  eruptions  also  have  triggered  tsunamis. 

2.  Tsunamis  are  principally  generated  by  undersea  earthquakes  of 
magnitudes  greater  than  6.5  on  the  Richter  scale.  The  typical  height  of 
a tsunami  in  the  deep  ocean  is  less  than  a foot  and  the  wave  period  is 

5 minutes  to  several  hours.  Tsunamis  travel  at  the  shallow-water  wave 
celerity  equal  to  the  square  root  of  acceleration  due  to  gravity  times 
water  depth  even  in  the  deepest  oceans  because  of  their  very  long  wave- 
lengths. This  speed  of  propagation  can  be  in  excess  of  500  mph*  in  the 
deep  ocean. 

3.  When  tsunami  waves  approach  a coastal  region  where  the  water 
depth  decreases  rapidly,  wave  refraction,  shoaling,  and  bay  or  harbor 
resonance  may  result  in  significantly  increased  wave  heights.  The  great 
period  and  wavelength  of  tsunami  waves  preclude  their  dissipating  energy 
as  a breaking  surf;  instead,  they  are  apt  to  appear  as  rapidly  rising 
water  levels  and  only  occasionally  as  bores. 

The  loss  of  life  and  destruction  of  property  due  to  tsunamis 
have  been  immense.  The  Great  Hoei  Tokaido-Nanhaido  tsunami  of  Japan 
killed  30,000  people  in  1707.  In  1868,  the  Great  Peru  tsunami  caused 

* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  and  metric  (SI)  units  to  U.  S.  customary 
units  is  presented  on  page  3. 


25,000  deaths  and  carried  the  frigate  U.S.S.  Waterlee  1,300  ft  inland. 

The  Great  Meiji  Sanriku  tsunami  of  1896  killed  27,122  persons  in  Japan 
and  washed  away  over  10,000  houses.  The  most  recent  major  tsunami  to 
affect  the  United  States,  the  196U  Alaskan  tsunami,  killed  107  people  in 
Alaska,  1*  in  Oregon,  and  11  in  Crescent  City,  California,  and  caused 
over  $100  million  in  damage  on  the  west  coast  of  North  America.1 

5.  The  Hawaiian  Islands,  a chain  of  eight  islands  as  shown  in 
Figure  1,  have  a history  of  destructive  tsunamis  generated  both  in  dis- 
tant areas  and  locally.  The  earliest  recording  of  a severe  tsunami  in 
the  Hawaiian  Islands  was  in  1837  when  a tsunami  from  Chile  reached  an 
elevation  of  20  ft  at  Hilo  and  killed  U6  people  in  the  Kau  District  of 
the  island  of  Hawaii  (Figure  l).  Prior  to  1837,  a number  of  severe 
tsunamis  undoubtedly  reached  the  islands,  but  unfortunately  no  detailed 

records  were  kept.  Since  1837,  there  have  been  16  tsunamis  that  have 

2 

caused  significant  damage. 

6.  The  most  destructive  tsunami  to  ever  hit  the  islands  in  terms 

of  loss  of  life  and  destruction  of  property  was  the  Great  Aleutian 
tsunami  of  19^6,  which  killed  173  people  and  produced  waves  over  55  ft 
in  elevation.  Hilo  incurred  $26  million  in  property  damage  attributable 
to  this  tsunami.  j \ j 

7-  The  i960  Chilean  tsunami  is  the  most  recent  distantly  gener- 
ated tsunami  that  produced  major  effects  in  the  Hawaiian  Islands.  Sixty- 
one  lives,  all  at  Hilo,  were  claimed  by  the  tsunami.  Damage  throughout 
the  State  was  estimated  to  be  $23.5  million  of  which  93  percent  occurred 
at  Hilo.  Other  major  damage  was  restricted  to  the  Kahului  area  of  the 
Island  of  Maui.  Inspection  of  the  damage  at  Hilo  revealed  much  evidence 
of  the  tremendous  forces  developed  by  the  waves.  Twenty-ton  boulders 
had  been  moved  hundreds  of  feet,  asphaltic  concrete  pavements  were 

peeled  from  their  subbase,  and  hundreds  of  automobiles  were  tossed 
3 

around  and  crushed. 

Tsunami  Inundation 

8.  Tsunami  inundation  is  the  flooding  during  part  of  a tsunami  j 
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event  of  land  that  is  normally  dry.  Inundation  is  used  in  this  report 
without  any  necessary  implication  of  severe  damage  or  great  water 
depths.  Dry  land  areas  flooded  during  a tsunami  are  in  the  region  of 
inundation. 

9.  Characteristics  of  tsunami  inundation  are  not  well  known 
since  there  have  been  few  surveys  of  both  tsunami  high-water  marks  near 
the  shoreline  and  detailed  measurements  of  inundation  limits.  Tsunamis 
in  the  Hawaiian  Islands  usually  appear  as  rapidly  rising  water  levels 
rather  than  bores.  Thus  runup  should  be  similar  to  elevation  at  the 
shoreline,  although  runup  may  be  greater  or  less  than  elevation  at  the 
shoreline  due  to  flow  divergence  or  convergence,  time-dependence  of  the 
wave  form,  frictional  effects,  and  topographic  variations. 

10.  The  similarity  between  runup  and  elevation  at  the  shoreline 

is  illustrated  in  the  detailed  surveys  that  have  been  performed.  For 

example,  Magoon  reports  flooding  to  about  the  20-ft  contour  above  mean 

lower  low  water  (mllw)  and  elevations  at  the  shoreline  of  about  20  ft 

(mllw)  for  the  196U  tsunami  at  Crescent  City,  California.  Wilson  and 

T^rum1  report  that  the  20-ft  (mllw)  runup  at  Valdez,  Alaska,  for  the 

196k  tsunami  checks  "well  for  consistency  with  water  level  measurements 

made  on  numerous  buildings  throughout  the  town."  Similar  comments  were 

made  by  Brown ^ in  reference  to  survey  measurements  of  30-ft  (mllw)  runup 

at  Seward,  Alaska,  for  the  1961*  tsunami.  Runup  and  elevation  at  the 

shoreline  w^re  similar  at  nine  locations  in  Japan  as  recorded  by  Nasu^ 

in  surveys  following  the  1933  Sanriku  tsunami.  This  tsunami  had  a short 

period  (12  minutes)  and  reached  an  elevation  as  great  as  28.7  metres  at 

one  survey  location.  Runup  and  elevation  at  the  shoreline  also  were 

similar  at  Hilo,  Hawaii,  for  the  i960  tsunami  (borelike  waves).  Dif- 

7 

ferences  are  apparent,  however,  at  locations  where  Eaton  et  al.  dem- 
onstrate that  flow  divergence  is  significant. 

11.  Tsunami  inundation  has  been  modeled  using  hydraulic 

8 9 

models.  ’ However,  hydraulic  models  must  be  built  large  enough  for 
scale  effects  not  to  be  significant.  Thus  it  is  quite  expensive  to 
build  hydraulic  models  for  simulations  of  tsunamis.  For  example,  al- 
though Hilo  Harbor  is  small  relative  to  the  entire  coastline  of  the 
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Hawaiian  Islands,  a hydraulic  model  of  this  harbor  that  was  used  for 
tsunami  studies  was  required  to  be  of  moderate  size  (right  triangular 
shape  with  sides  measuring  63  and  93  ft)  even  with  a small  horizontal 

g 

scale  of  1:600  and  a three  to  one  distortion. 

12.  There  have  been  several  analytical  and  numerical  studies  of 
one-dimensional  tsunami  runup. Comparisons  were  not  made  in  these 
studies  of  predicted  runup  and  actual  inundation  during  historical 
tsunamis.  The  one-dimensional  nature  of  the  solutions  and  also  the 
borelike  waves  considered  in  the  studies  make  these  studies  inappro- 
priate for  most  practical  applications.  A recent  study  by  Bretschneider 

ll+ 

and  Wybro  employs  a one-dimensional  analytic  solution  to  determine 
tsunami  runup  for  both  bores  and  nonbores.  Frictional  effects  and 
linear  ground  slopes  are  included  in  the  solution.  Comparisons  are 
presented  of  calculated  and  measureu  runup  of  actual  historical  tsunamis. 
The  main  drawbacks  of  this  approach  are  the  one-dimensional  and  time- 
independent  properties  of  the  solution. 

13.  A one-dimensional  empirical  method  based  upon  trial-and-error 
comparisons  of  calculated  (assuming  various  initial  elevations  and  rates 
of  elevation  decay  during  propagation  over  land)  and  measured  runup  dur- 
ing historical  tsunamis  in  Hawaii  has  been  presented  by  Cox.^  Cox 
sought  criteria  by  which  the  long-term  limits  of  potential  inundation 
might  be  estimated  with  necessary  safety  margins.  Prior  to  this  work 
the  Civil  Defense  of  the  State  of  Hawaii  warned  people  to  evacuate  areas 
where  the  land  elevation  was  less  than  50  ft  above  sea  level.  In  some 
areas  such  as  Honolulu,  such  an  evacuation  area  would  have  been  much 
larger  than  could  have  been  evacuated  or  would  have  needed  to  have  been 
evacuated.  Cox  found  that  an  assumed  initial  total  energy  head  of  50  ft 
(at  a location  where  the  water  depth  was  10  ft  below  msl)  and  a one  per- 
cent decline  with  distance  traversed  (l-ft  decay  for  every  100  ft  of 
travel)  would  provide  an  estimate  of  the  potential  inundation  limits 
that  appeared  reasonable  when  compared  with  historical  records.  Some- 
what different  criteria  were  applied  at  certain  other  locations. 

ll+.  Two-dimensional  numerical  models  have  not  been  developed  to 
simulate  tsunami  inundation.  The  pioneering  work  in  numerical  models 
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for  flooding  problems  was  a study  of  storm  surges  in  Galveston  Bay, 

Texas,  by  Reid  and  Bodine.1^  An  extension  of  this  work  was  performed  by 
17 

Masch  et  al.  and  applied  to  tidal  hydrodynamics  problems. 

Purpose  of  Study 

15.  The  purpose  of  this  study  was  to  develop  a method  to  deter- 
mine two-dimensional  and  time-dependent  tsunami  (nonbore)  inundation 

based  upon  a numerical  model  developed  recently  at  the  WES  for  numerical 

l8 

simulations  of  tidal  hydrodynamics.  Such  a method  is  needed  by  the 

U.  S.  Army  Engineer  Division,  Pacific  Ocean  (POD),  to  delineate  inunda- 

near  the  shoreline  pre- 
19 

WES  report . Inundation 
flood  hazard  evaluations 
rate  calculations. 


tion  limits  of  tsunamis  based  upon  elevations 
dieted  for  the  Hawaiian  Islands  in  a previous 
limits  are  required  by  POD  for  use  in  tsunami 
for  floodplain  management  and  flood  insurance 
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PART  II:  NUMERICAL  MODEL 


Equations  of  Motion 


16.  The  hydrodynamic  equations  governing  tsunami  propagation  in 
the  nearshore  region  and  inundation  are  derived  from  the  classical 

20 

Navier-Stokes  equations  expressed  in  a Cartesian  coordinate  system. 

Since  tsunamis  have  lengths  much  greater  than  the  water  depths  over 
which  they  propagate,  fluid  motions  are  approximately  two-dimer sional. 
Thus  vertical  fluid  velocities  and  accelerations  are  small  in  comparison 
with  horizontal  velocities  and  accelerations,  respectively.  If  the 
fluid  is  assumed  to  be  homogeneous  and  incompressible  and  the  flow  is 
integrated  from  the  sea  bottom  to  the  free  water  surface,  the  usual  long 
wave  equations  are  obtained: 


JW  + U.3U  V3U  ajii  gU 

8t  d 3x  d 3y  sa3x  2 2 


3V  + U 3V  + V 3V 
3t  d 3x  d 3y 


ff,3n  + _gv_ 

6 * c2d2 


(u" 


V2) 


1/2 


= o 


V2) 


1/2 


= 0 


(1) 


(2) 


In  + 9U  + 3V 
3t  3x  3y 


0 


(3) 


where 

U = vertically  integrated  transport  per  unit  width  in 
x direction* 

t = time 

d = n - h , total  water  depth 

h = still-water  depth 

x,y  = Cartesian  coordinates 

V = vertically  integrated  transport  per  unit  width  in 
y direction 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de- 
fined in  the  Notation  (Appendix  A). 
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g = acceleration  due  to  gravity 
n = water-surface  elevation 
C = Chezy  frictional  coefficient 

1/2  1/2 

IT.  The  terms  (U2  + V2)  and  (U2  + V2)  are 

C d (Td 

bottom  friction  terms  derived  from  experiments  and  theoretical  consider- 

21  22  23 

ations  for  one-dimensional  flow  in  a river.  Previous  efforts  ’ in 
numerical  modeling  of  long  waves  have  shown  that  these  terms  adequately 
represent  the  damping  of  long  waves  due  to  the  shear  at  the  sea  bottom. 
10.  The  Chezy  frictional  coefficient  is  given  by  the  expression 


C = 


1.U86  .1/6 

d 


(M 


where  nm  is  Manning's  n . Manning's  n is  empirically  determined 
and  depends  upon  characteristics  of  the  land  over  which  water  is  flowing. 

19.  The  advective  terms  in  Equations  1 and  2 (second  and  third 
terms)  are  neglected  in  this  study.  These  terms  often  cause  difficul- 
ties in  numerical  modeling  since  they  generate  higher  harmonics  of  waves 
and  intensify  vorticity  through  stretching  of  vortex  lines.  The  genera- 
tion of  higher  harmonics  may  create  problems  since  the  cascading  of 
energy  from  longer  to  shorter  waves  will  result  in  energy  accumulation 
in  waves  with  a wavelength  of  twice  the  grid  spacing  of  the  numerical 
grid.  This  accumulation  may  result  in  numerical  instability.  Difficul- 
ties in  modeling  the  generation  and  dissipation  of  vorticity  also  can  re- 
sult in  eddy  formations  that  may  produce  numerical  instability. 

20.  The  advective  terms  are  not  important  in  determining  the  ex- 
tent of  tsunami  inundation.  These  terms  were  found  to  have  a negligible 
effect  on  the  extent  of  tsunami  inundation  in  tests  (with  and  without 
advective  terms)  of  the  numerical  model  presented  later  in  this  report. 
Accurate  numerical  simulations  of  hurricane  surge  and  tidal  elevations 

have  been  made  in  the  past,  solving  equations  that  neglect  the  advective 

. 2k 

terms. 


! 
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Variable  Grid 


21.  Tsunami  inundation  is  difficult  to  model  since  tsunamis  have 
great  wavelengths.  Wavelengths  even  in  shallow  water  are  much  greater 
than  extents  of  inland  flooding  or  length  scales  of  bathymetric  and 
topographic  variations  near  the  coastline.  The  input  boundary  of  a 
numerical  grid  must  be  at  least  one-half  a wavelength  away  from  the 
coastline,  since  the  form  of  the  waves  reflected  from  the  coast  is  not 
known  a priori  at  the  input  boundary.  Thus  the  forcing  function  at  the 
input  boundary  cannot  be  varied  to  account  for  both  the  incident  and 
reflected  wave.  This  makes  it  necessary  to  move  the  input  boundary  far 

enough  away  from  the  coast  that  waves  are  not  reflected  at  this  boundary. 

1 8 

22.  The  numerical  model  employed  in  this  study  uses  a smoothly 
varying  grid  that  allows  cells  to  be  small  in  the  inundation  region  and 
large  in  the  ocean.  A piecewise  reversible  transformation  (analogous  to 
that  used  in  Reference  25)  is  used  independently  in  the  x-  and  y-direc- 
tions  to  map  the  variable  grid  into  a uniform  grid  used  in  the  computa- 
tional space.  The  transformation  has  the  form 

x = a1  + b^a^Cl  (5) 

y = a2  + *>2a2c2  (6) 

where  a^,b^,c^,a2,b2,c2  are  arbitrary  constants  and  and  a2  are 

coordinates  in  the  computational  space.  This  transformation  allows  all 
derivatives  to  be  centered  in  the  computational  space.  Many  stability 
problems  commonly  occurring  in  variable  grid  schemes  are  eliminated 
using  this  transformation  since  the  real  space  grid  is  smoothly  varying 
with  the  variation  and  its  first  derivative  being  continuous. 

23.  The  equations  of  motion  (Equations  1-3)  have  the  following 
form  in  the  computational  space: 
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O SURFACE  ELEVATION  <tj)  — FLOW/UNIT  WIDTH  IN 

WATER  DEPTH  (d)  x-DIRECTION  (U) 

FRICTIONAL  COEFFICIENT  (C)  1 y-OIRECTION '(V)H 


Figure  2.  Space-staggered  scheme 


where  d#  = n - h and  At  is  the  time-step  for  completing  two  cycles. 
A single  bar  over  a variable  represents  a two-point  average  and  a double 
bar  a four -point  average.  The  subscripts  m and  n correspond  to 
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spatial  locations  and  the  superscript  k to  time  levels.  If  terms 
are  collected  that  are  to  he  computed  along  line  n at  time  level 
(k  + l/2)At  , Equations  12  and  13  can  he  expressed  as  follows: 


I . I 


-a  nk+1^2  + £ nk+1/2  + k+1/2 

m n,m  m+l/2  n, m+l/2  am+l  ^n.m+l  ~ m+l/2 


,*+1/2  k+1/2  TTk+l/2  _ . 

m-1/2  n,m-l/2  ^n,m  m+l/2  n, m+l/2  m 


where 


At 

= n 


Aa2  I n+l/2,m  ~ n-l/2,m 


R = ,*-1/2  + a+  I Kd»  ) k-1/2  k-l/2 

m+l/2  n,m+l/2  I 2^y^  Aal  |r'n>m+l  ” nn,m 


® n, m+l/2 


2C2(d*)2 


The  coefficients  are  defined  as  follows: 


Uk~1/2  V 

n, m+l/2 j 


Vl/2  ~ am+l/2  = 2^1jAa1 


a = a = gAtd* 

m m+1  2^ij  Aa^ 


r r 


, , gAt 
am+l/2  " 1 .-2, 


2C‘:(d»)‘ 


(■ 


nk-1/2  f + y2 

n,m+l/2 i V 


1/2 


(20) 


k+l/2 

25.  Assuming  that  n is  a given  water  level  at  the  lower 

k+l/2  ^ 

boundary  and  U L+1/2  is  a 8iven  flow  at  the  uPPer  boundary,  the  matrix 
form  of  the  equations  for  line  n can  be  written  as  follows: 


Si+1/2  ***4+1  0 ° 

~aM+l/2  1 ***4+3/2  0 

II  I I I 

- n . (2D 

0 -aM+l  aM+3/2  aM+2  * ’ * 0 


• . -a. 


UM+l/2 

k+l/2 

BM+l/2 

nM+l 

\l+ 1 

UM+3/2 

= 

BM+3/2 

nL 

L-l/2 


26.  The  above  system  of  equations  is  tridiagonal  and  can  be 
solved  with  a minimum  number  of  operations  by  performing  the  process  of 
elimination  using  a set  of  recursive  formulas.  Equation  1^  is  expressed 
as  follows: 


R = R + « nk+l/2 

m+1/2  Bm+l/2  m nM 


Substitution  of  Equation  21  into  the  Equation  15  yields 


m+1/2 


(- B.  C/2  * * C/2  * V3/2  C3/2  - Vi  w 


where 


T^+1/2  = -P  Uk+1/2  + 0 
m+1  m+1  m+3/2  Vl 


P = m+3/2 

m+1  1 + a , R 

m+1/2  m 


0 = Am+1  + am+l/2  Sm 

m+1  1 + a . , ,_R 

m+1/2  m 


Again,  the  flow  rate  can  be  expressed  as  a function  of  the  next  water 
level  using  Equation  lit. 

^+3/2  = "Rm+1  nm+2  + Sm+1  ^ 


where 


am+3/2  + am+lRm+l 


_ Bm+3/2  + am+lSn+l 


&m+3/2  + am+lRm+l 


In  general,  the  recursion  formulas  can  be  written  as  follows: 


IT 


(32) 


where 


where 


T1  = 1 + a , R , 
m-1/ 2 m-1 


P _ m+l/2 

m - T1 


T2  = a . /o  + a P 
m+l/2  m m 


Am  + am-l/2  Sm-1 


27.  Fractional  subscripts  are  not  permitted  in  the  FORTRAN 
computer  code.  Therefore,  in  the  computer  code  a different  index  system 
is  adopted.  Water  levels,  depths,  frictional  coefficients,  and  flows 
appearing  in  the  cell  shown  in  Figure  2 have  the  same  coordinate  index 
(N,M)  . Using  this  notation,  the  coefficients  given  in  Equations  26, 

27,  and  29-30  are  expressed  as  follows: 


28.  The  solution  Equations  21  and  31  can  be  expressed  as  follows: 


= -rm  . njj+jj/2  + S 

N,M-1  M-l  N,M  M-l 


k+1/2  p yk+1/2  . 

nN,M  M N,M  74 


29.  The  recursion  coefficients  of  Equations  3^-37  can  be  computed 
in  succession  between  boundaries  on  line  N . Various  approximations  of 
these  coefficients  are  required  depending  upon  the  types  of  boundary 
conditions  that  are  applied.  The  solution  Equations  and  h6  can  be 
solved  for  all  surface  elevations  and  flows  in  descending  order  once 
the  set  of  coefficients  has  been  calculated. 

30.  The  second  cycle  computes  n and  V implicitly,  advancing 
the  time  from  (k  + l/2)At  to  (k  + l)At  . The  recursion  formulas  for 
this  cycle  are  very  similar  to  those  given  for  the  first  cycle.  These 
formulas  for  the  second  cycle  (applied  along  a line  m)  can  be 
expressed  as  follows: 


Vn+l/2,m  = Vn+l/2,m  + At  ( " 2^  [^j  ^ {W 


k+1  . k k (.  SVn+l/2,m  + 

- n + n - n /+  , 

n,m  n+l,m  n,ml  / 


^•yk+l  + ^ 

V n+l/2,m  n+l/2,m/ 

(2  i c42 


^n+l/2,m 
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irf- 


where 

C _ -k+1/2  - 

d = n - h 

31.  The  recursion  formulas  for  the  second  cycle  can  be  written  in 
the  same  notation  as  Equations  3i*-46. 


B 


N 


(55) 


a ■ 1/2(v^M  * <M/2  - Vl,»  - Vm)  (56) 


u = i /u(u^+1/2  + uk+1^2  + uk+1^2  + uk+1/,2\  (57) 

U X/^\UN+1,M  N+l.M-l  UN,M-1  UN,M  J V?" 


32.  The  solution  Equations  21  and  31  can  be  expressed  as  follows: 


r)k+'''  = -P  Vk+1  + Q. 
nN,M  N N,M  TJ 


(58) 


Vk+1  - _P  nk+1  + Q 

N-1,M  “ rl-l  nN,M  N-l 


(59) 


Boundary  Conditions 

33.  There  are  four  types  of  boundaries  used  in  the  numerical 
model.  The  first  is  a seaward  boundary  that  terminates  the  computational 
grid.  Water  levels  are  prescribed  as  a function  of  location  and  time  at 
this  boundary.  Thus,  water  levels  may  vary  along  the  input  boundary. 
Water  level  time-histories  can  either  be  arbitrary  with  elevations  sup- 
plied in  tabular  form  or  sinusoidal  with  the  amplitude  and  period 
prescribed. 

3*+.  A second  boundary  used  is  located  at  the  land-water  interface. 
The  no-flow  condition  normal  to  this  boundary  is  the  usual  boundary  con- 
dition employed  in  the  numerical  computations.  Thus  U = 0 or  V = 0 
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is  prescribed  at  the  appropriate  grid  cell  face.  The  model  also  handles 
the  flooding  and  drying  of  land.  Inundation  is  simulated  by  making  the 
location  of  land-water  boundary  a function  of  local  water  elevations. 

The  possibility  of  inundation  is  established  by  continual  monitoring  of 
water  levels  in  adjacent  grid  cells.  Initial  movement  of  water  onto 
previously  dry  grid  cells  is  controlled  by  a broad-crested  weir 
formula.1^  Once  the  water  level  on  a previously  dry  cell  exceeds  some 
small  prescribed  value,  the  boundary  face  is  treated  as  being  open  and 
computations  for  n , U , and  V are  the  same  as  for  any  ocean  cell. 
The  drying  of  cells  follows  the  inverse  process. 

35.  Exposed,  submerged,  and  overtopping  barriers  defined  along 
grid  cell  faces  form  a third  type  of  boundary.  Exposed  barriers  have 
no-flow  conditions  prescribed  across  appropriate  grid  cell  faces.  Flow 
across  a submerged  barrier  is  controlled  by  use  of  a time-dependent 
frictional  coefficient.  Overtopping  barriers  are  barriers  that  are 
exposed  (with  no-flow  conditions)  during  times  of  a simulation  and  sub- 
merged at  other  times.  The  overtopping  of  a barrier  is  controlled  by 
using  a broad-crested  weir  formula  to  specify  proper  flow  rates  across 
the  barrier.  Water  is  transferred  from  the  high  to  low  side  using  this 
formula.  Once  the  barrier  is  submerged  (or  exposed)  the  procedures  used 
for  submerged  (or  exposed)  barriers  are  applied. 

36.  A no-flow  condition  normal  to  the  boundary  is  used  for  the 
lateral  ocean  boundaries  of  the  computational  grid.  The  waves  reflected 
from  the  shoreline  have  only  small  propagation  vector  components  in  a 
direction  perpendicular  to  these  boundaries.  There  is  some  reflection 
from  these  boundaries;  however,  the  reflected  waves  still  move  away  from 
the  shoreline  until  they  reach  the  input  boundary.  Computations  are 
stopped  before  reflections  from  these  lateral  boundaries  are  a problem. 
The  lateral  boundaries  may  introduce  small  errors  in  the  immediate 
vicinity  of  the  intersection  of  these  boundaries  with  the  land-water 
interface. 
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PART  III:  VERIFICATION 


Introduction 


37-  There  have  been  fe«  extensive  surveys  of  tsunami  high-water 

marks  near  the  shoreline  and  detailed  measurements  of  inundation  levels 

in  addition  to  recordings  of  tsunami  wave  forms.  The  only  two  tsunamis 

in  the  United  States  for  which  such  data  are  available  are  the  i960 

Chilean  tsunami  recorded  at  Hilo,  Hawaii,  and  the  1961*  Alaskan  tsunami 

recorded  at  Crescent  City,  California.  However,  the  incident  wave  form 

for  the  i960  tsunami  at  Hilo  was  quite  complex.  The  major  waves  of  this 
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tsunami  at  Hilo  were  bores.  Wiegel  indicates  that  Mach  reflection  of 
a wave  may  have  occurred  on  the  cliffs  of  the  west  side  of  Hilo  Bay  with 
the  Mach-stem  wave  moving  toward  the  Hilo  area.  This  Mach-stem  wave  may 
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have  superposed  on  the  diffracted  incident  wave  at  Hilo.  Tsunami  ele- 
vations recorded  near  the  shoreline  for  this  tsunami  vary  quite  rapidly 
along  the  shoreline.  The  i960  tsunami  at  Hilo  was  not  used  for  verifica- 
tion purposes  in  this  study,  since  the  incident  waves  for  this  tsunami 
were  bores;  the  complex  Mach-stem  effect  may  have  developed;  and  waves 
traveling  in  different  directions  may  have  superposed  with  unknown 
phases.  The  lack  of  knowledge  concerning  the  incident  and  reflected 
waves  would  make  simulations  of  this  tsunami  difficult  even  if  the  model 
could  handle  bores. 

38.  The  196U  tsunami  at  Crescent  City,  California  (Figure  3),  was 
used  to  verify  the  numerical  model  discussed  in  PART  II.  The  simula- 
tion of  this  tsunami  is  an  excellent  verification  test  since  the  tsunami 
was  quite  large  at  Crescent  City  and  the  region  has  many  complex  fea- 
tures. For  example,  the  Crescent  City  Harbor  is  protected  by  break- 
waters, some  of  which  were  overtopped  and  others  which  were  not.  In  the 
region  there  is  a developed  city  area,  mud  flats,  and  an  extensive 
riverine  floodplain.  Inland  flooding  was  extensive  in  the  floodplain 
area  and  extended  as  much  as  a mile  inland.  Sand  dunes  and  elevated 
roads  played  a prominent  role  in  limiting  flooding  in  certain  areas. 
Verification  of  the  model  for  this  complex  region  makes  the  model  quite 
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Figure  3.  Location  map  for  Crescent  City, 
California 


general  so  that  it  can  be  applied  to  any  region. 

1* 

39-  Magoon  reports  surveys  of  the  Crescent  City  region  made  by 
the  Corps  of  Engineers  following  the  196**  tsunami . The  maximum  limits  of 
inundation  and  high-water  marks  were  recorded  in  the  surveys.  Wiegel2^ 
presents  characteristics  of  the  wave  form  and  details  of  the  interaction 
of  the  tsunami  with  the  Crescent  City  region. 

Numerical  Grid  and  Boundary  Conditions 


(mllw).  Beyond  this  contour,  the  ocean  is  a constant  30-ft  depth  (mllw;. 
This  constant  depth  region  allows  the  input  boundary  to  be  moved  suffi- 
ciently distant  (at  least  one-half  the  total  wave  form  of  interest)  from 
the  shoreline  that  waves  do  not  reflect  from  the  input  boundary.  The 
variable  nature  of  the  grid  allows  the  input  boundary  to  be  placed  far 
from  the  region  of  interest  without  creating  so  many  grid  cells  that 
computational  time  becomes  excessive.  A uniform  grid  with  cells  the 
size  of  the  smallest  cells  in  Figure  k would  require  a computational 
time  of  approximately  25  times  greater  than  that  required  for  the  vari- 
able grid  shown  in  Figure  4. 

1+1.  The  input  to  the  numerical  model  is  a wave  crest  that,  when 
propagated  to  shore,  reproduces  as  accurately  as  possible  the  historical 
maximum  elevation  at  the  tide-gage  location  near  the  shoreline  in 
Crescent  City  Harbor.  The  tide  gage  was  destroyed  during  the  1961+ 
tsunami',  however,  Wilson  and  T^rum"*-  inferred  the  form  of  the  wave  from 
the  tide  gage  recording  prior  to  destruction  in  addition  to  other  evi- 
dence. There  is  some  disagreement  among  investigators  as  to  the  maxi- 
mum elevation  of  the  tsunami  at  the  tide  gage.  Estimates  range  from  18 
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to  20.7  ft  above  mllw.  However,  according  to  Wiegel,  the  20.7-ft  ele- 
vation was  recorded  on  land  at  a distance  of  approximately  1000  ft  from 

the  tide  gage  and  it  included  a 1-  to  2-ft  contribution  from  local  runup 
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on  a concrete  wall.  An  earlier  study  by  Keulegan  et  al.  considered 
all  available  information  and  concluded  that  the  maximum  elevation  at 
the  tide  gage  was  19  ft  above  mllw.  A wave  crest  with  an  elevation  of 
19  ft  above  mllw  at  the  tide-gage  location  was  used  in  the  numerical 
simulation. 

1+2.  There  are  several  barriers  in  the  Crescent  City  region. 

27 

Wiegel  reports  that  there  was  little  flow  over  the  outer  breakwater 
(Figure  5)>  Therefore,  this  breakwater  was  represented  as  a non- 
overtopping barrier.  However,  the  breakwaters  attached  to  Whaler  Island 
(Figure  5)  were  overtopped  during  a part  of  the  tsunami  and  thus  are 
represented  as  dynamic  overtopping  barriers.  Similarly,  Redwood  Highway 
(Figure  5)  and  beach  dunes  along  the  crest  were  overtopped  during  the 
1961+  tsunami.  Water  flowed  over  these  barriers  in  a manner  similar  to 
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PACIFIC  OCEAN 

Figure  5.  Crescent  City  Harbor  area 

li 

water  flowing  over  a broad  weir  and  the  runup  was  reduced.  These 
barriers  also  are  represented  as  dynamic  overtopping  barriers  with  flow 
rates  governed  by  weir  equations. 

1+3.  Frictional  effects  in  the  numerical  model  are  governed  by  the 

Chezy  frictional  coefficient  given  in  Equation  U.  A spatially  varying 

Manning's  n must  be  selected  at  the  center  of  each  grid  cell  shown  in 

Figure  h.  Manning's  n is  a function  of  the  roughness  of  the  ground 

and  the  vegetation.  Values  of  Manning's  n can  be  found  in  numerous 

references.  ’ In  this  study,  values  of  Manning's  n suggested  by 
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Bretschneider  and  Wybro  for  various  coastal  terrain  conditions  were 
used.  Manning's  n values  were  selected  for  very  general  categories  of 
terrain,  since  detailed  knowledge  of  vegetation  and  land  roughness  are 
usually  not  known  for  an  area.  For  the  Crescent  City  region,  the  ocean 
bottom,  mud  flats,  and  beaches  were  assigned  a Manning's  n equal  to 
0.021+.  Developed  areas  were  assigned  a value  of  0.035  and  the  riverine 
floodplain  area  and  other  heavily  vegetated  areas  a value  of  0.055.  No 


attempt  was  made  to  force  agreement  of  numerical  calculations  and  histor- 
ical recordings  of  elevations  by  varying  local  values  of  Manning's  n . 

UU.  Ocean  depths  and  land  elevations  were  taken  from  the  Sister 
Rocks,  California,  quadrangle  map  published  by  the  U.  S.  Geological 

Survey.  Elevations  in  the  city  of  Crescent  City  were  taken  from  the 
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detailed  surveys  of  the  Corps  of  Engineers  as  reported  by  Magoon.  Land 
elevations  are  not  known  accurately  in  the  large  floodplain  area  (Elk 
Creek)  near  Crescent  City  so  elevations  had  to  be  estimated  in  this  area. 

Results 

^5.  The  Corps  of  Engineers  made  surveys  of  the  tsunami  inundation 
at  Crescent  City  in  January,  April,  and  May  of  1965  as  reported  by 
Magoon.  Figure  6 shows  the  location  of  11  high-water  marks  recorded  by 
the  surveys.  The  tide  gage  was  located  at  the  end  of  Citizen's  Dock 


Figure  6.  High-water  mark  locations  for  196U  tsunami  (adapted  from 

Magoon,  1965) 
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(Figure  5)  near  high-water  mark  number  307.  The  following  tabulation 
presents  the  measured  elevations  at  the  high-water  marks  and  elevations 
calculated  by  the  numerical  model.  The  measured  and  calculated  eleva- 
tions are  within  1 ft  except  at  locations  307,  2,  3,  and  5-  However, 
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Wiegel  reports  that  the  measured  elevation  at  location  307  was  greater 
than  the  actual  wave  elevation  by  1 to  2 ft  due  to  local  runup  on  a con- 
crete wall.  Thus  the  actual  wave  elevation  and  the  calculated  elevation 
are  in  good  agreement  at  location  307.  The  measured  elevation  at  loca- 
tion 4 is  lower  than  the  calculated  elevation  since  the  measured  high- 
water  mark  was  inside  a lumber  building.  Incomplete  filling  of  the 
lumber  building  by  water  would  result  in  lower  elevations  within  the 
building  than  those  outside  the  building.  Location  number  5 is  located 


High-Water 

Mark  Number 
(Reference  4) 

Measured 
Elevation 
mllw,  ft 

Calculated 
Elevation 
mllw,  ft 

Difference 

ft 

302 

18.35 

17.95 

+0.40 

305 

18.74 

19-39 

-0.65 

307* 

20.70 

19.45 

+1.25 

309 

19.84 

19.79 

+0.05 

312 

19.41 

20.23 

-0.82 

316 

16.29 

17.09 

-0.80 

1 

15.90 

15.43 

+0.47 

2 

17.80 

19.24 

-1.44 

3 

19.30 

18.89 

+0.41 

16.50 

18.94 

-2.44 

5+ 

20.50 

19.13 

1.37 

Note:  Mean  error  (all  locations)  = 0.2  ft;  root  mean  square  error  (all 
locations)  = 1.1  ft;  root  mean  square  error  (locations  307,  4,  5 
eliminated)  = 0.7  ft. 

* Measured  elevation  1-2  ft  greater  than  actual  wave  elevation  as 
a result  of  runup  on  concrete  wall  (Reference  27). 

**  Measured  elevation  inside  lumber  building  (Reference  4).  Eleva- 
tion greater  outside  building. 

t Coastal  location  almost  2 miles  south  of  tide  gage.  Incident 
wave  amplitude  likely  to  be  somewhat  different  here  than  at  tide 
gage  location. 


almost  2 miles  south  of  the  tide  gage.  It  is  likely  that  the  wave  ap- 
proaching this  shoreline  had  a greater  elevation  than  the  wave  approach- 
ing the  harbor  area.  The  elevation  calculated  at  this  location  could  be 
increased  by  varying  the  incident  wave  form  along  the  input  boundary. 
There  is  no  apparent  explanation  for  the  disagreement  at  location  2, 
especially  since  there  is  good  agreement  at  nearby  location  1. 

46.  The  agreement  between  the  measured  high-water  marks  and  cal- 
culated elevations  is  quite  good.  The  root-mean-square  error  is  0.7  ft 
(elevations  at  locations  307,  4,  and  5 not  included  since  differences 
attributable  to  factors  other  than  the  numerical  model).  Since  the  ob- 
served elevation  at  the  tide  gage  used  to  choose  an  incident  wave  height 
is  accurate  to  within  no  more  than  0.5  to  1 ft,  the  agreement  between 
measured  and  calculated  elevations  is  good. 

47.  Figure  7 shows  contour  lines  of  the  tsunami  elevation  above 
ground  level  within  the  developed  area  of  Crescent  City.  The  topography 
and  tsunami  elevations  are  known  in  much  more  detail  in  this  developed 


Figure  7.  Inundation  lines,  1964  tsunami  (adapted  from  Magoon,  1965) 
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area  than  any  other  area.  The  dotted  parts  of  the  measured  elevation 
dashed  lines  denote  uncertain  elevations.  The  contour  lines  for  the 
calculated  elevations  were  determined  by  linear  interpolation.  For 
example,  if  the  elevation  at  the  center  of  one  cell  was  9 ft  above  land 
and  the  elevation  in  sin  adjacent  cell  was  7 ft,  the  contour  line  was 
drawn  half  way  between  the  two  cell  centers.  The  agreement  between  the 
measured  and  calculated  elevations  shown  in  Figure  7 is  good. 

48.  Figure  8 shows  inundation  contours  for  the  complete  region 
surveyed  by  the  Corps  of  Engineers.  The  most  accurate  contours  are  in 
the  developed  region  of  Crescent  City.  Contours  in  floodplain  areas  are 
of  unknown  accuracy,  since  the  surveys  were  performed  approximately  one 
year  after  the  tsunami.  It  may  have  been  difficult  to  accurately  deter- 
mine  inundation  limits  in  the  floodplain  areas  since  there  were  no 
structures  in  these  areas,  debris  lines  may  have  been  obliterated  or 
difficult  to  differentiate  from  river-swept  debris,  and  vegetation  may 


Figure  8.  Flooding  for  1964  tsunami  (adopted  from  Magoon,  1965) 
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have  recovered  from  salt  damage  as  a result  of  the  freshwater  flow  of 
the  river  system  (winter  rains  occurred  between  the  time  of  the  tsunami 
and  the  surveys ) . 

1+9-  The  measured  and  calculated  inundation  lines  in  the  city  area 
(Figure  8)  are  in  excellent  agreement.  The  calculated  inundation  lines 
are  farther  inland  than  the  measured  lines  in  the  floodplain  areas . In 
addition  to  the  question  of  the  accuracy  of  the  measured  inundation  in 
the  floodplain  areas,  the  topography  of  the  floodplains  is  not  known 
accurately.  The  elevations  of  the  floodplains  are  only  known  to  be 
greater  than  mean  sea  level  (msl)  and  less  than  10  ft  above  msl.  Ele- 
vations throughout  this  area  had  to  be  estimated  in  the  topographic 
input  to  the  numerical  model.  Elevations  calculated  at  the  extent  of 
flooding  in  these  floodplain  areas  were  small  (often  a fraction  of  a 
foot).  Land  elevations  could  have  been  increased  by  small  increments  to 
have  prevented  flooding  beyond  the  measured  limits  and  thus  forced 
agreement  between  measured  and  calculated  limits.  However,  this  was  not 
done,  since  the  accuracy  of  the  survey  measurements  in  these  areas  was 
questionable  and  the  topography  of  floodplain  areas  is  rarely  accurate. 
Neither  friction  coefficients  nor  land  elevations  were  varied  in  this 
verification  to  force  agreement  with  measured  elevations,  since  in  an 
application  of  the  numerical  model  to  an  arbitrary  location  these  param- 
eters would  not  be  accurately  known.  The  good  agreements  shown  in  Fig- 
ure 7 and  the  tabulation  on  page  30  indicate  that  small  inundation 
levels  in  distant  floodplain  areas  do  not  have  significant  influence  on 
elevations  in  other  areas. 
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PART  IV:  APPLICATION  TO  HAUULA-PUNALUU , HAWAII 

Introduction 

50.  The  tsunami  inundation  numerical  model  was  developed  to  allow 
the  POD  to  determine  inundation  limits  when  tsunami  elevations  near  the 
shoreline  are  known.  Elevations  near  the  shoreline  produced  by  tsunamis 
with  various  return  periods  are  known  for  the  entire  coastline  of  the 
Hawaiian  Islands  (except  the  coast  of  the  uninhabited  U.  S.  Navy  target 
island  of  Kahoolawe)  from  Reference  19-  As  an  example  of  the  use  of  the 
numerical  model,  the  model  was  applied  to  the  Hauula-Punaluu  region  of 
the  Island  of  Oahu  (Figure  9)  to  determine  flooding  produced  by  50-  and 
100-year  tsunamis. 

51.  Topographic  data  and  descriptions  of  vegetation  covering  the 
land  for  the  Hauula-Punaluu  region  were  taken  from  a flood  hazard  map 
of  the  region  published  by  the  State  of  Hawaii  in  cooperation  with  the 
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POD.  This  map  has  a scale  of  1:4800.  Bathymetric  data  for  the  region 
were  taken  from  the  Hauula  and  Kahana  quadrangle  maps  of  the  U.  S.  Geo- 
logical Survey  (scale  of  1:24,000). 

Numerical  Grid  and  Boundary  Conditions 

52.  Figure  10  shows  the  numerical  grid  used  in  the  simulation  of 
tsunami  inundation  in  the  Hauula-Punaluu  region.  The  grid  is  oriented 
such  that  the  incident  wave  approaches  the  Hauula-Punaluu  region  from  a 
direction  perpendicular  to  the  60-ft  (mllw)  bathymetric  contour.  It  was 
assumed  that  refraction  effects  beyond  the  60-ft  contour  bend  the  propa- 
gation path  of  a tsunami  such  that  it  approaches  the  nearshore  area 
moving  in  a direction  approximately  perpendicular  to  this  contour.  The 
ocean  bathymetry  is  reproduced  out  to  the  60-ft  contour  (mllw)  and  be- 
yond this  contour  the  ocean  is  assumed  to  be  a constant  60-ft  depth 
(mllw) . 

53.  The  sizes  of  the  grid  cells  in  the  region  that  is  above  msl 
were  selected  to  allow  a good  representation  of  the  land  topography. 

Cells  are  longer  in  the  direction  approximately  parallel  to  shoreline 
than  in  the  direction  perpendicular  to  it  since  land  elevations  vary 
more  gradually  in  a direction  parallel  to  the  shoreline  (land  contours 
are  approximately  parallel  to  the  shoreline).  Cell  sizes  in  the  near- 
shore ocean  region  were  selected  for  the  same  reason. 

54.  The  input  to  the  numerical  model  is  a wave  that  when  propa- 
gated to  shore  reproduces  desired  elevations  200  ft  inland  from  shore- 
line at  sites  in  Hauula  and  Punaluu.  A distance  of  200  ft  inland  from 
the  shoreline  was  selected  since  the  elevations  presented  in  Refer- 
ence 19  were  based  upon  historical  data  recorded  approximately  200  ft 
inland  (on  the  average).  Locations  of  these  two  sites  are  presented  in 
Reference  33.  Flooding  was  determined  for  elevations  that  are  equaled 
or  exceeded  once  in  100  years  and  once  in  50  years  at  the  two  sites. 
Reference  19  predicts  100-  and  50-year  elevations  of  12.6  and  9-7  ft, 
respectively,  at  Hauula  and  13.0  and  10.4  ft,  respectively,  at  Punaluu. 

55.  A trial-and-error  method  is  used  to  produce  desired  elevations 
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at  the  Hauula  and  Punaluu  sites.  Different  incident  waves  are  tested 
until  maximum  elevations  attain  desired  values  at  locations  that  are 
200  ft  inland  from  these  two  sites.  The  numerical  model  allows  the  inci- 
dent wave  form  to  vary  along  the  input  boundary.  However,  the  great  dis- 
tance between  the  input  boundary  and  the  shoreline  that  was  selected  in 
order  to  eliminate  the  problem  of  reflected  waves  allows  enough  time 
during  propagation  for  energy  to  move  along  the  wave  front  and  reduce 
the  variability.  To  alleviate  this  problem  it  was  necessary  to  add  a 
thin  nonpermeable  barrier  between  the  two  sites  that  extended  from  the 
input  boundary  along  the  length  of  the  constant  depth  region.  This 
method  allows  an  input  boundary  to  be  located  in  effect  near  the  shore- 
line without  the  problem  of  reflected  waves  striking  the  true  input 
boundary.  Since  the  barrier  is  parallel  to  the  direction  of  propagation 
of  the  tsunami,  it  does  not  create  a problem  by  reflecting  waves  back 
toward  the  shoreline. 

5 6.  There  were  seven  different  Manning's  n values  used  to  repre- 
sent frictional  effects  in  the  Hauula-Punaluu  region.  The  ocean  bottom 
and  beaches  were  assigned  a Manning's  n equal  to  0.02H.  Developed 
areas  were  assigned  a Manning's  n of  0.35,  bush-  and  tree-covered 
areas  or  rough  floodplain  areas  a value  of  0.055,  cultivated  areas  or 
floodplain  areas  without  bushes  and  trees  a value  of  0.0^5,  marsh  areas 
a value  of  0.0^0,  reefs  a value  of  0.150,  and  noncultivated  areas  with- 
out bushes  and  trees  a value  of  0.030. 

Results 
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57-  Figures  11-20  show  contours  of  50-  and  100-year  tsunami  inun- 
dations calculated  using  the  numerical  model  for  the  Hauula-Punaluu  re- 
gion extending  from  Kaipapau  Point  to  Makalii  Point.  The  figures  have 

32 

the  same  scale  as  the  map  referenced  in  the  introduction  to  PART  IV  of 
this  report.  Thus  the  figures  may  be  overlaid  on  this  map  to  delineate 
detailed  extents  of  inundation. 

58.  Inundation  contours  were  determined  first  by  delineating 
cells  in  the  land  area  that  were  flooded  during  a part  of  the  tsunami 
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Figure  11.  Tsunami  inundation  in  che  vicinity  of  fCaipapau  Point 
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Figure  12.  Tsunami  inundation  in  the  vicinity  of  Hauula  Beach  Park 
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Figure  13.  Tsunami  inundation  in  the  vicinity  of  Hauula 
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Figure  lU.  Tsunami  inundation  west  of  Kalaipaloa  Point 
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Figure  15.  Tsunami  inundation  east  of  Kalaipaloa  Point 
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Figure  l6.  Tsunami  inundation  south  of  Kaluanui  Stream 
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Figure  17.  Tsunami  inundation  in  the. vicinity  of  Papaakuko 
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Figure  18.  Tsunami  inundation  north  of  Punaluu 
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Figure  19-  Tsunami  inundation  in  the  vicinity  of  Punaluu 
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Figure  20.  Tsunami  inundation  in  the  vicinity  of  Makalii  Point 


and  also  bordered  cells  that  did  not  flood.  These  cells  that  did  flood 
would  have  depths  of  water  upon  them  that  were  not  great  enough  to  flood 
onto  adjacent  cells  that  remained  dry  during  the  tsunami.  However, 
there  would  be  enough  water  to  flood  to  some  intermediate  elevation  be- 
tween the  centers  of  the  flooded  and  dry  cells.  For  example,  if  the 
elevation  at  the  center  of  one  cell  was  1*+  ft  and  it  was  covered  with 
1 ft  of  water  and  the  elevation  at  the  center  of  the  adjacent  dry  cell 
was  17  ft,  inundation  in  the  prototype  would  be  to  approximately  the 
15-ft  contour  between  the  cell  centers.  Thus  for  the  length  of  this 
cell  the  level  of  inundation  would  be  the  15-ft  contour.  Since  spatial 
variations  of  individual  contours  are  typically  less  than  cell  sizes , 
the  inundation  curves  shown  in  Figures  11-20  often  have  spatial  varia- 
tions smaller  than  typical  cell  sizes. 

59.  The  map  of  Reference  32  plots  inundation  contours  for  a 
100-year  tsunami.  A tsunami  elevation  of  15  ft  at  the  shoreline  was 
selected  for  this  map  by  increasing  recorded  elevations  in  the  region  for 
the  19^6  tsunami  by  25  percent.  It  was  then  assumed  that  the  energy  of 
the  tsunami  would  dissipate  at  a rate  that  would  result  in  a declining 
maximum  water-surface  profile  of  1 ft  for  every  100  ft  inland.  However, 
such  a decline  in  water-surface  profile  is  much  more  rapid  than  has  been 
observed.  For  example,  Hatori  reports  a decline  an  order  of  magnitude 
smaller  based  upon  three  different  tsunamis.  The  one  percent  decay  used 
by  Cox^  applies  to  the  total  energy  head  and  not  the  water-surface 
level.  The  100-year  tsunami  inundation  contour  predicted  in  this  study 
is  farther  inland  than  the  100-year  contour  presented  in  Reference  32  at 
all  locations  except  for  a small  section  in  Figure  12. 

60.  Figure  21  shows  contours  of  tsunami  elevations  above  the  land 
surface  for  the  100-year  tsunami  in  the  vicinity  of  Kalaipaloa  Point. 

This  figure  illustrates  the  complicated  two-dimensional  pattern  of  flood- 
ing. The  Poki  Wai  Place  area  was  flooded  to  a depth  a little  less  than 

6 ft  since  it  is  a residential  area  that  is  built  up  to  a higher  eleva- 
tion than  neighboring  land.  Poki  Wai  Place  is  surrounded  by  a canal 
filled  with  water. 

61.  The  numerical  model  predicts  some  interesting  elevation 
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Figure  21.  Tsunami  elevations  above  the  land  surface 
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patterns.  For  example,  at  the  shoreline  there  is  an  elevation  "hump," 
that  is,  elevations  are  larger  at  the  shoreline  than  just  inland  or 
seaward  of  the  shoreline.  This  increase  in  elevation  may  be  due  to  the 
reflection  of  a wave  from  the  shoreline  and  the  subsequent  superposition 
of  incident  and  reflected  waves,  shoaling  of  the  incident  wave,  or  con- 
version of  kinetic  energy  to  potential  energy  as  the  increased  friction 
slows  the  incident  waves.  There  also  are  many  instances  where  land  at 
the  same  elevation  experiences  different  levels  of  inundation.  Some- 
times land  floods  at  one  location  but  not  at  another  location  that  has 
the  same  elevation.  Such  inundation  differences  may  be  attributable  to 
flow  divergence  and  convergence,  two-dimensional  variations  in  roughness, 
reflections  by  barriers,  and  topographic  variations.  For  example,  land 
near  a riverine  floodplain  may  not  flood  as  a result  of  flow  passing  the 
land  to  flood  the  lower  lying  floodplain  area. 


PART  V:  CONCLUSIONS 


62.  The  finite  difference  numerical  model  presented  in  this 
report  accurately  predicts  tsunami  high-water  marks  and  inundation  for 
nonbore  tsunamis.  The  model  simulates  two-dimensional,  time-dependent, 
and  dissipative  effects.  A conformal  coordinate  transformation  allows 
the  grid  to  be  variable  so  that  small  grid  cells  can  be  concentrated  in 
areas  where  high  resolution  is  required.  The  grid  cells  can  be  large  in 
the  ocean  region  where  high  resolution  is  not  necessary,  since  bathy- 
metric variations  are  gradual  and  tsunami  wavelengths  are  long.  The 
variable  grid  makes  it  feasible  to  calculate  tsunami  inundation  over  a 
large  area  since  it  is  not  necessary  to  use  very  small  grid  cells  over 
the  entire  region. 

63.  The  tsunami  inundation  numerical  model  was  developed  to  be 

quite  general  so  that  it  could  be  used  to  predict  tsunami  flooding  at 

arbitrary  locations.  The  model  allows  nonpermeable  breakwaters, 

permeable  breakwaters,  dynamic  overtopping  breakwaters,  elevated  roads, 

sand  dunes,  and  other  overtopping  barriers.  Frictional  dissipation, 

land  topography,  and  the  ocean  bathymetry  can  vary  in  any  arbitrary  two- 

dimensional  fashion.  Since  tsunamis  do  not  usually  appear  as  bores,  the 

model  is  applicable  to  most  tsunami  inundation  problems.  There  also  is 
15  7 

evidence  * ’ that  inundation  produced  by  bores  is  similar  to  inundation 
produced  by  nonbores. 

6U.  The  numerical  model  presented  in  this  report  is  an  engineer- 
ing tool  that  the  POD  can  use  to  determine  tsunami  inundation  throughout 

the  Hawaiian  Islands.  The  model  has  been  documented  in  a separate  WES 
35 

report . 
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APPENDIX  A:  NOTATION 

al 

Constant  coefficient 

a2 

Constant  coefficient 

a 

m 

Coefficient  defined  at  location  m 

Am 

Variable  defined  at  location  m 

bl 

Constant  coefficient 

b2 

Constant  coefficient 

B 

m 

Variable  defined  at  location  m 

ci 

Constant  coefficient 

C2 

Constant  coefficient 

C 

Chezy  frictional  coefficient 

d 

Total  water  depth,  ft 

d* 

d 

__ 

n - h , averaged  total  water  depth,  ft 
-k+l/2  - 

n - h , averaged  total  water  depth,  ft 

» 

D1 

Variable 

g 

2 

Acceleration  due  to  gravity,  32.2  ft/sec 

h 

Still-water  depth,  ft 

\ 

k 

Integral  associated  with  time  levels 

L 

Integral  grid  index 

m 

Integral  grid  index 

M 

Integral  grid  index 

n 

Integral  grid  index 

n 

m 

Manning’s  n 

Pm 

Variable  defined  at  grid  index  m 

Sn 

Variable  defined  at  grid  index  m 

R 

m 

Variable  defined  at  grid  index  m 

S 

m 

Variable  defined  at  grid  index  m 

t 

Time,  sec 

At 

Time-step  for  completing  two  cycles , sec 

T1 

Variable 

T2 

Variable 

|r 

[ 

Vertically  integrated  transport  per  unit  width  in  x-direction, 
ft^/sec 

Al 

1 

- 

l 

• 

V Vertically  integrated  transport  per  unit  width  in  y-direction, 
ft^/sec 

x Cartesian  coordinate,  ft 

y Cartesian  coordinate,  ft 

Coordinate  in  computational  space 

atp  Coordinate  in  computational  space 

n Water-surface  elevation,  ft 

Variable  defining  stretching  of  computational  grid 
1^2  Variable  defining  stretching  of  computational  grid 
t Time-step  for  completing  one  cycle,  sec 


Two-point  time  average  when  over  variable 
Four-point  time  average  when  over  variable 


’ 


V 

1 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
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